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Abstract 

We show that a tetragonal lattice of weakly interacting cavities with uniaxial electromagnetic 
response is the photonic counterpart of topological crystalline insulators, a new topological phase 
of atomic band insulators. Namely, the frequency band structure stemming from the interaction 
of resonant modes of the individual cavities exhibits an omnidirectional band gap within which 
gapless surface states emerge for finite slabs of the lattice. Due to the equivalence of a topological 
crystalline insulator with its photonic-crystal analog, the frequency band structure of the latter can 
be characterized by a topological invariant. Such a topological photonic crystal can be realized 
in the microwave regime as a three-dimensional lattice of dielectric particles embedded within a 
continuous network of thin metallic wires. 

PACS numbers: 42.70.Qs, 73.20.-r, 73.43.-f 
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I. INTRODUCTION 



The frequency band structure of artificial periodic dielectrics formally known as photonic 
crystals is the electromagnetic (EM) counterpart of the electronic band structure in ordi- 
nary atomic solids. Recently, a new analogy between electron and photon states in periodic 
structures has been proposed by Raghu and Haldane,- namely the one-way chiral edge states 
in two-dimensional (2D) photonic-crystal slabs which are similar to the corresponding edge 
states in the quantum Hall effect.- The photonic chiral edge states are a result of time- 
reversal (TR) symmetry breaking which comes about with the inclusion of gyroelectric/ 
gyromagnetic material components; these states are robust to disorder and structural im- 
perfections as long as the corresponding topological invariant (Chern number in this case) 
remains constant. 

In certain atomic solids, TR symmetry breaking is not prerequisite for the appearance 
of topological electron states as it is the case in the quantum Hall effect. Namely, when 
spin-orbit interactions are included in a TR symmetric graphene sheet, a bulk excitation 
gap and spin-filtered edge states emerge^ without the presence of an external magnetic field, 
a phenomenon which is known in literature as quantum spin Hall effect. Its generalization 
to three-dimensional (3D) atomic solids lead to a new class of solids, namely, topological 
insulators. 4 The latter possess a spin-orbit-induced energy gap and gapless surface states 
exhibiting insulating behavior in bulk and metallic behavior at their surfaces. 

Apart from topological insulators where the spin-orbit band structure with TR symmetry 
defines the topological class of the corresponding electron states, other topological phases 
have been proposed such as topological superconductors (band structure with particle-hole 
symmetry),- magnetic insulators (band structure with magnetic translation symmetry) ,- 
and, very recently, topological crystalline insulators.- In the latter case the band structure 
respects TR symmetry as well as a certain point-group symmetry leading to bulk energy 
gap and gapless surface states. 

In this work, we propose a photonic analog of a topological crystalline insulator. Our 
model photonic system is a 3D crystal of weakly interacting resonators respecting TR sym- 
metry and the point-symmetry group associated with a given crystal surface. As a result, 
the system possesses an omnidirectional band gap within which gapless surface states of 
the EM field are supported. It is shown that the corresponding photonic band structure is 
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equivalent to the energy band structure of an atomic topological crystalline insulator and, 
as such, the corresponding states are topological states of the EM field classified by a Z 2 
topological invariant. 

The frequency band structure of photonic crystals whose (periodically repeated) con- 
stituent scattering elements interact weakly with each other can be calculated by a means 
which is similar to the tight-binding method employed for atomic insulators and semicon- 
ductors. Photonic bands amenable to a tight-binding-like description are e.g., the bands 
stemming from the whispering-gallery modes of a lattice of high-index scattered the defect 
bands of a sublattice of point defects, within a photonic crystal with an absolute band gap,- 
the plasmonic bands of a lattice of metallic spheres^ or of a lattice of dielectric cavities within 
a metallic host.— In the latter case, the frequency band structure stems from the weak in- 
teraction of the surface plasmons of each individual cavity 1 ^ wherein light propagates within 
the crystal volume by a hopping mechanism. Such type of lattice constitutes the photonic 
analog of a topological crystalline insulator presented in this work whose frequency band 
structure will be revealed based on a photonic tight-binding treatment within the framework 
of the coupled-dipole method.— The latter is an exact means of solving Maxwell's equations 
in the presence of nonmagnetic scatterers. 

II. TIGHT-BINDING DESCRIPTION OF DIELECTRIC CAVITIES IN A PLAS- 
MONIC HOST 

We consider a lattice of dielectric cavities within a lossless metallic host. The i-th cavity 
is represented by a dipole of moment P« = (Pi- X , Pi- y , P^ z ) which stems from an incident 
electric field E mc and the field which is scattered by all the other cavities of the lattice. This 
way the dipole moments of all the cavities are coupled to each other and to the external 
field leading to the coupled-dipole equation 



Gu' (u) is the electric part of the free-space Green's tensor and cti(u) is the 3x3 polarizability 
tensor of the z-th cavity. Eq. (pQ) is a 3iV x 3N linear system of equations where N is 
the number of cavities of the system. We assume that the cavities exhibit a uniaxial EM 
response, i.e., the corresponding polarizability tensor is diagonal with and 




(1) 
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ol z = ot±. For strong anisotropy, the cavity resonances within the xy-p\ane and along the 
z-axis can be spectrally distinct; thus, around the region of e.g., the cavity resonance Un 
within the xy-plane, a± <C aii (see appendix). In this case, one can separate the EM response 
within the xy-plane from that along the z-axis and Eq. ([T]) becomes a 2N x 2N system of 
equations, 

P i = a l |(w)[^G«/(a;)P P ]. (2) 

where we have set E mc = since we are seeking the eigenmodes of the system of cavities. 
Also, now, Pj = (/',;,. / , , :?; ). 

For a particle/cavity of electric permittivity e\\ embedded within a material host of per- 
mittivity €h, the polarizability an is given by the Clausius-Mossotti formula 

3V e \\ - e h /ox 

an = 3 

4tt e|| + 2e h K J 

where V is the volume of the particle/ cavity. For a lossless plasmonic (metallic) host in 
which case the electric permittivity can be taken as Drude-type, i.e., eh — 1 —Up/u 2 (where 
Up is the bulk plasma frequency), the polarizability ay exhibits a pole at = u p ^j2/ (ejj + 2) 
(surface plasmon resonance). By making a Laurent expansion of ocn around ouu and keeping 
the leading term, we may write 

«ll = — = 7j ( 4 ) 
where F = (w||/2)(e|| — e/i) / + 2). For sufficiently high value of the permittivity of the 
dielectric cavity, i.e., en > 10, the electric field of the surface plasmon is much localized at 
the surface of the cavity. As a result, in a periodic lattice of cavities, the interaction of 
neighboring surface plasmons is very weak leading to much narrow frequency bands. By 
treating such a lattice in a tight binding-like framework, we may assume that the Green's 
tensor Gu>(oj) does not vary much with frequency and therefore, Gh'(uj) ~ G&'(u>[|). In this 
case, Eq. (j2J) becomes an eigenvalue problem 

J2Gii'(^\\)Pi'=^i (5) 

where 



Gni\u)\\ 



C(<?||N'|)I 2 + J{q\\\r iV \) 



x u'Vn' 



x u>yu> Vji> 



(6) 



with ra> = Ti — r»/, q\\ = y/thUJ\\/c and I2 is the 2x2 unit matrix. The form of functions 
C(qu \ru'\), J(qu \ru>\) generally depends on the type of medium hosting the cavities (isotropic, 
gyrotropic, bi-anisotropic, etc).— ^ The Green's tensor of Eq. describes the electric in- 
teractions between two point dipoles Pj and each of which corresponds to a single cavity. 
The first term of Gu' describes an interaction which does not depend on the orientation of 
the two dipoles whilst the second one is orientation dependent. 
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FIG. 1: (Color online) (a) Tetragonal crystal with two cavities within the unit cell, (b) The bulk 
Brillouin zone and (c) the surface Brillouin zone corresponding to the (001) surface. 

For an infinitely periodic system, i.e., a crystal of cavities, we assume the Bloch ansatz 
for the polarization field, i.e., 

Pi = P n/3 = exp(zk ■ R n )P 0/ 3 (7) 

The cavity index i becomes composite, i = nf3, where n enumerates the unit cell and (3 
the positions of inequivalent cavities in the unit cell. Also, R n denotes the lattice vectors 
and k = (k x , k y , k z ) is the Bloch wavevector. By substituting Eq. ([7]) into Eq. ([5]) we finally 
obtain 

^G^(w||,k)P ^ = nP ^ (8) 

0' 

where 

G^'(^||,k) = ^exp[ik • (R„ - Rn'))G n fc n ,p,(u\\). (9) 
Solution of Eq. ([8]) provides the frequency band structure of a periodic system of cavities. 



III. TOPOLOGICAL FREQUENCY BANDS 

Since Eq. ([8]) is equivalent to a Hamiltonian eigenvalue problem, we adopt the crystal 
structure of Ref. Q. Namely, a tetragonal lattice with a unit cell consisting of two same 
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FIG. 2: Frequency band structure for tetragonal lattice of resonant cavities within a plasmonic 
host (see Fig. [I]) corresponding to the Green's tensor of Eq. (|lip with = — s 2 = 1.2, s 2 = — s 2 = 
0.5,5; = 2.5,4 = °- 5 ^ = 2 - 

cavities at inequivalent positions A and S [see Fig. QJa)] along the c-axis. In this case, the 
index /3 in Eq. (jSJ) assumes the values /3 = A,B for each sublattice (layer) of the crystal. The 
above lattice is characterized by the C4 point-symmetry group. In order to preserve the C4 
symmetry^ in the Green's tensor matrix of Eq. (jSJ) we assume that the interaction between 
two cavities within the same layer (either A or B) depends on the relative orientation of 
the point dipole in each cavity whilst the interaction between cavities belonging to adjacent 
layers is orientation independent. Also, we take into account interactions up to second 
neighbors in both inter- and intra- layer interactions. Taking the above into account, the 
lattice Green's tensor assumes the form 




(10) 



where 




G AB {k) = [s[ + 2s' 2 (cos(k x a) + cos(k y a)) + s' z exp(ik z a)]l2- 



(11) 
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FIG. 3: Frequency band structure for a finite slab ABAB ■ ■ ■ ABB of the crystal of Fig. [TJ made 
from 80 bilayers. 

The lattice Green's tensor of Eq. ( TTTT) is completely equivalent to the lattice Hamiltonian 
of Ref. sf , s 2 , s' 2 i s z in Eq. (Hip generally depend on gy, the lattice constant a and 
the interlayer distance c but hereafter will be used as independent parameters. Namely we 
choose sf = -sf = 1.2, s 2 = —sf = 0.5, = 2.5, s' 2 = 0.5, s z = 2. In Fig. [2] we show the 
(normalized) frequency band structure corresponding to Eq. (|TT|) along the symmetry lines 
of the Brillouin zone shown in Fig. QJb). It is evident that an omnidirectional frequency 
band gap exists around Q = which is prerequisite for the emergence of surface states. 
In order to inquire the occurrence of surface states we find the eigenvalues of the Green's 
tensor of Eq. (fTTj) in a form appropriate for a slab geometry. The emergence of surface 
states depends critically on the surface termination of the finite slab, i.e., for different slab 
terminations different surface-state dispersions occur (if occur at all). Namely, we assume a 
finite slab parallel to the (001) surface (characterized by the C4 symmetry group) consisting 
of 80 alternating AB layers except the last bilayer which is BB, i.e., the layer sequence is 
ABAB ■ ■ ■ ABB. The corresponding frequency band structure along the symmetry lines of 
the surface Brillouin zone of the (001) surface [see Fig. [U^c)] is shown in Fig. |3j It is evident 
that there exist gapless surface states within the band gap exhibiting a quadratic degeneracy 
at the M-point. In this case, the corresponding doublet of surface states can be described 
by an effective theory^ similarly to the doublet states at a point of linear degeneracy (Dirac 
point).— 

We note that the equivalence of the Green's tensor G with the atomic Hamiltonian of 
Ref. Q as well as the form of the time-reversal T and (geometric) C*4-rotation U operators 
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for the EM problem-* which are the same as for spinless electrons, allows to describe the 
photonic band structure with the Z 2 topological invariant v$ 



\'A) 



\VTM ( 1 \ VAZ 



(12) 



where for real eigenvectors of G(k) we have 7 



4) klk2 = Pf [ W (k 2 )]/Pf ^(ki)] 



(13) 



and 

Wmn(kj) = (u m (k i )|E/'|u n (k i )). (14) 

Pf stands for the Pfaffian of a skew-symmetric matrix, i.e., Pf [w] 2 = det(io). Due to the 
double degeneracy of the band structure at the four special momenta points T, M, A, Z the 
frequency bands come in doublets. Since frequency eigenvectors with different eigenfrequen- 
cies are orthogonal, all the inter-pair elements of the w-matrix are zero and the latter is 
written as: 

/ 



w(ki 







^ 2 (k t ) 
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y w N {k i ) J 

where w J (kj) are anti-symmetric SU(2) matrices,— i.e., u> J (kj) = A\ or A2, where 



(15) 



Ai 



1 
-1 



A, 



-1 

1 



(16) 



In this case, Pf[u>(kj)] = w\ 2 w\ 2 • ■ "^12 = Therefore, (— l) klk2 = ±1 and vq = 1 which 
ensures the presence of gapless surface states. 

We note that the above analysis relies on the assumption of real frequency bands. The 
presence of losses in the constituent materials renders the frequency bands complex, i.e., the 
Bloch wavevector possesses both a real and an imaginary part. However, even in this case, 
one can still speak of real frequency bands if the imaginary part of the Bloch wavevector is 
at least hundred times smaller than the corresponding real part. This a common criterion 
used in calculations of the complex frequency band structure by on-shell electromagnetic 
solvers such as the layer- multiple scattering method^ or the transfer-matrix method. 2 - 
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FIG. 4: (Color online) A possible realization of a photonic structure with gapless surface states: 
dielectric particles of square cross section, joined together with cylindrical coupling elements and 
embedded within a 3D network of metallic wires (artificial plasma). 

IV. BLUEPRINT FOR A PHOTONIC TOPOLOGICAL INSULATOR 

A possible realization of the photonic analogue of topological insulator in the laboratory is 
depicted in Fig. |H Since our model system requires dielectric cavities within a homogeneous 
plasma, a lattice of nano-cavities formed within a homogeneous Drude-type metal, e.g., a 
noble metal (Au, Ag, Cu), would be the obvious answer.— However, the plasmon bands 
are extremely lossy due to the intrinsic absorption of noble metals in the visible regime. 
A solution to this would be the use of an artificial plasmonic medium operating in the 
microwave regime where metals are perfect conductors and losses are minimal. Artificial 
plasma can be created by a 3D network of thin metallic wires of a few tens of /zm in diameter 
and spaced by a few mm.— A lattice of dielectric particles within an artificial plasma can 
be modelled with the presented tight-binding Green's tensor. Since the interaction among 
first and second neighbors within the same bilayer (A or B) should depend on the dipole 
orientation (in order to preserve the C4 symmetry), the dielectric particles in each layer are 
connected with cylindrical waveguiding elements (different in each layer A or B - see Fig. HI). 
In contrast, between two successive bilayers there are no such elements since interactions 
between dipoles belonging to different layers should be independent of the dipole orientations 
(s orbital-like). Another advantage of realizing the photonic analog in the microwave regime 
is the absence of nonlinearities in the EM response of the constituent materials since photon- 
photon interactions may destroy the quadratic degeneracy of the surface bands in analogy 
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with fermionic systems.— 

Finally, we must stress that a photonic topological crystalline insulator can be also re- 
alized with purely dielectric materials if the host medium surrounding the cavities is not a 
plasmonic medium but a photonic crystal with an absolute band gap: the cavities would 
be point defects within the otherwise periodic photonic crystal and the tight-binding de- 
scription would be still appropriate.- 1 ^ In this case, Maxwell's equations lack of any kind of 
characteristic length and the proposed analog would be realized in any length scale. 



V. CONCLUSIONS 



In conclusion, a 3D lattice of weakly interacting cavities respecting TR symmetry and 
a certain point-group symmetry constitutes a photonic analog of a topological crystalline 
insulator by demonstrating a spectrum of gapless surface states. A possible experimental 
realization would be a 3D lattice of dielectric particles within a continuous network of thin 
metallic wires with a plasma frequency in the GHz regime. 

This work has been supported by the European Community's Seventh Framework Pro- 
gramme (FP7/2007-2013) under Grant Agreement No. 228455-NANOGOLD (Self-organized 
nanomaterials for tailored optical and electrical properties). 



Appendix 

The ^-component of the polarizability need not be zero but can assume finite values as 
far as the frequency band stemming from the surface-plasmon resonance corresponding to 
the z-component is spectrally distinct (no overlap) with the bands stemming from the xy- 
components of the polarizability. In this case, one can treat separately the two frequency 
bands (doublet) stemming from the xy-resonance from the band coming from the z-resonance 
(singlet). The above requirements can be quantified as follows, 

c, « ^± (A.i) 

where Cm is the width of the xy- frequency bands (in dimensionless frequency units). Cn 
is obtained from the first term of the right-hand side of Eq. (|6]) of the paper. To a first 
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approximation, it is given by~ 



^ exp(- g „q,|) 

Q||a]| 



where gy = y|eftja;||/c. Therefore, the condition (1A.1I) is written as 



giiaii ca 



where a\\ is the lattice constant in the xy-plane. Given that u\\ = u p y/2/{e\\ + 2), Eq. (1A.3|) 
becomes 



exp(- v /|e h |w||a||/c) ^ / 2(e ± - ey) 



|e h |w||0||/c ~ V ( e ll+ 2 )(^ + 2 ) 
From the above equation it is evident that for a given value of the dielectric anisotropy ej_— en, 
one can always find a suitably large lattice constant a\\ such that Eq. (1A.4I) is fulfilled. The 
latter allows the easy engineering of the photonic analog of a crystalline topological insulator 
since there is practically no restriction on the choice of the (uniaxial) material the cavities 
are made from. It can also be easily understood that if Eq. ([6]) holds, the same equation is 
true for the width C± of the singlet frequency band (resulting from the z- resonance) . 

Numerical example. Suppose that the cavities are made from a nematic liquid crystal 
which is a uniaxial material. Typical values of the permittivity tensor e are e.g., en = 1.5, 
e_i_ = 1.8. In this case, u\\ ~ 0.75u p and = 1 — uj^/uj ~ —0.777. By choosing a large 
lattice constant, i.e., a\\ = Ac/u p , Eq. ( 1A.4I) is clearly fulfilled 



exp (— \/\eh\u)\\a\\/c) / 2(e i — en) , , , 

Vl 1 0.026866 < 0.2213 w J- ^7 — ( A - 5 ) 



e h \uj\\a\\/c - — y + 2 )( e± + 2) 
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